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Abstract. We here discuss the process of opinion formation in an open community where agents are made 
to interact and consequently update their beliefs. New actors (birth) are assumed to replace individuals 
that abandon the community (deaths). This dynamics is simulated in the framework of a simplified model 
that accounts for mutual affinity between agents. A rich phenomenology is presented and discussed with 
reference to the original (closed group) setting. Numerical findings are supported by analytical calculations. 

PACS. 87.23. Ge Dynamics of social systems - 05.45.-a Nonlinear dynamics and nonlinear dynamical 
systems 

1 Introduction drive the process of opinion making [I] . This large pro- 

duction has not been always accompanied by an adequate 
Opinion dynamics modeling represents a challenging field benchmarking effort to the relevant psychological litera- 
where ideas from statistical physics and non-linear sci- ture; and doubt s are consequently being cast on the in- 
ence can be possibly applied to understand the emergence terpretative ability of the proposed mathematical formu- 
of collective behaviors, like consensus or polarization in lations . Despite these intrinsic limitations, several models 
social groups. Several toy models have been proposed in 
the past to reproduce the key elements that supposedly 
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of interacting agents display a rich and intriguing dynam- Indeed, the outcome of an hypothetic binary interaction 

ics which deserves to be fully unraveled. relies on the difference of opinions, previously postulated, 

^ . . , . , , ,,.„,. , but also on the quality of the mutual relationships. The 

Opinion dynamics models can be classified m two large 

,, , , . . j i • affinity is dynamically coupled to the opinion, and, in this 

groups. On the one hand, opinions are represented as dis- J J 

z . ,., \ , ,, , , respect, it introduces a memory bias into the system: affin- 

crete (spm-like) variables where the system behaves sim- r ' J J 

• nl , -i , i m n ,i, ,. i ■ ity between agents increases when their opinion tends to 

narly to spin glasses models [2j. On the other, each m- J ° ^ 

dividual bears a continuous opinion which span a pre- converge. 

, rm T , , , , The aforementioned model is here modified to accom- 

assigned range 3J. In both approaches, a closed system 

j . , -, , , c modate for a death/born like process. In this formulation, 

is generally assumed, meaning that the same pool of ac- ' i ' 

, . j • i-i , , . mi • M agents are randomly eliminated from the system, ev- 

tors is made to interact during the evolution. Ihis can ° J ' 

... , , ,, , . , n . , ery T time steps. When an agent exits from the commu- 

be interpreted by assuming that the inherent dynamical 

, , . . % , r nity (virtually, dies), he is immediately replaced by a new 

timescales (e.g. opinion convergence time) are much faster J v J ' " j r j 

, j j 1 1 / . element, whose opinion and affinity with respect to the 

than those associated to the processes (e.g. migration, ' * j r 

i-i / j ,,s r j-c .. C ii group are randomly assigned. As we shall see, the pertur- 

birth/death) responsible for a modification of the group ° r jo i ±- 

, ,, rr i i • i • ii r bation here prescribed alters dramatically the behavior of 
composition, these latter enects having being therefore so 

r lijoi t ■ ■ -i the system, with reference to the ideal close-system config- 

lar neglected, buch an implicit assumption is certainly cor- J ' jo 

, , „ r uration. To understand such modifications via combined 

rect when the debate is bound to a small community of 

. ,. . , , , . ., , , . numerical and analytical tools, constitutes the object of 

individuals, thus making it possible to eventually achieve 

, c , a the investigations here reported, 

a rapid convergence towards the final conhguration. Oon- ° 

i , • i, i -i The paper is structured as follows. We first introduce 

versely, it might prove inaccurate when applied to a large 

, , r ■ ,i i the model, then we present the obtained analytical and 

ensemble of interacting agents, as the process becomes r ■' 

., , , , j, i u jj-u numerical results and, finally, we sum up and draw our 

considerably slower and external perturbations need to be 

accounted for. Given the above, it is therefore of interest corLCms i ons - 

to elucidate the open system setting, where the population 

is periodically renewed. ^ The model 

To this end, we refer to the model presented in [3], In the following, we will shortly review the model previ- 

where the role of affinity among individuals is introduced ously introduced in [4] and present the additional features 

as an additional ingredient. This novel quantity measures that are here under inspection. The interested reader can 

the degree of inter-personal intimacy and sharing, an ef- thus refer to the original paper [1] for a more detailed 

feet of paramount importance in real social system [5] . account on the model characteristics. 
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Consider a population of N agents and assume that at When the two agents i and j are selected on the basis 
time t they bear a scalar opinion 0\ <G [0, 1]. We also intro- of the recipe prescribed above, they interact and update 



duce the NxN time dependent matrix a , whose elements their characteristics according to the following scheme 
a' - belong to the interval [0, 1]. The quantities a\ ■ specify 



the affinity of individual i vs. j, at time t: Larger values \ (3) 

alf 1 = a\, + a*, [1 - af,]r 2 {AO l3 ) , 

of alj are associated to more trustable relationships. K 

Both the affinity matrix and the agents opinions are where the functions A and T 2 respectively read: 

randomly initialized time t = 0. At each time step t, two , t , 1 r „. t , n 

J\ (a*,-) = - [©(^ - a c ) + 1] (4) 



2 

and 



agents, say i and J, are selected according to the following 

extraction rule: first the agent i is randomly selected, with 

a uniform probability. Then, the agent j which is closer to ^ 2 = I ~ AO c ) (5) 

i in term of the social metric DV. is selected for interaction. and the symbol stands for the Heaviside step-functionl]. 

The quantity D% results from the linear superposition of More spec ifi ca ii y] p x is or 1 while P 2 is -1 or 1, de- 

the so-called social distance, d%j , and a stochastic contri- pending on the value of their respective arguments. In the 

bution rjj, namely: above expressions a c and AO c are constant parameters. 

n i7 _ rt , / n \ /i\ Notice that, for a c — * 0, the opinion and affinity are for- 

mally decoupled in ([3]), and the former evolves following 

the Deffuant et al. scheme [3] with /x = 0.5. 

In [1] , a complete analysis of the qualitative behavior 

of the model as a function of the involved parameters is 

reported. Asymptotic clusters of opinion are formed, each 

d\j = AOlj(l — alj) j = 1, TV j i , (2) agglomeration being different in size and centered around 

, , , , distinct opinion values. Individuals sharing the same be- 

with AO\ s = \0\-0)\. 



Here 77^ (0, cr) represents a normally distributed noise, 
with mean zero and variance cr, the latter being named 
social temperature. The social distance is instead defined 
as: 



1 The evolution of the quantities Oj (t) and ctij (t) is straight- 
The smaller the value of dy the closer the agent j 

forwardly obtained by switching the labels i and j in the equa- 

to i, both in term of affinity and opinion. The additive 

tions. 

noise ?7j(0,(t) acts therefore on a fictitious ID manifold, 2 In thg switchers A and p 2 are smoot h functions con- 
which is introduced to define the pseudo-particle (agent) s tmcted from the hyperbolic tangent. We shall here limit the 
interaction on the basis of a nearest neighbors selection discussion to considering the Heaviside approximation, which 
mechanism and, in this respect, set the degree of mixing is recovered by formally sending /3i,2 to infinity in Eqs. (3) and 
in the community. (4) of [4]. 



4 



Timoteo Carletti et al.: Opinion Dynamics in an Open Community 





Fig. 1. Left panel: Typical evolution of the opinion versus time, i.e. number of iterations and the asymptotic distribution of 
agents' opinion in a closed community. Right panel: Final affinity matrix, represented by the underlying network: node i is linked 
to node j if ay > a = 0.8. Here a = 4 • 10 -4 , AO c = 0.5, a c — 0.5 and p — 0, i.e. no agent can leave the group. Initial opinion 
are (random) uniformly distributed within the interval [0, 1], while a^j is initialized with uniform (random) values between 
an 0.5. 



lieves are also characterized by a large affinity scores, as 
it is exemplified in Figure [TJ 

More quantitatively, the system is shown to undergo a 
continuous phase transition: above a critical value of the 
control parameter (erac) -1 / 2 the system fragments into 
several opinion clusters, otherwise convergence to a sin- 
gle group is numerically shown to occur. The interested 
reader can refer to [4] for a speculative interpretation of 
this finding in term of its psychological implications. We 
shall here simply notice that a significant degree of mixing 
(large social temperature a) brings the system towards the 
single-cluster final configuration. 

Starting from this setting, we introduce the birth/death 
process, which in turn amounts to place the system in con- 
tact with an external reservoir. The perturbation here hy- 
pothesized is periodic and leaves the total number of agent 
unchanged. Every T time steps (i.e. encounter events) M 



agents, randomly selected, are forced to abandon the sys- 
tem (death). Every removed individual is instantaneously 
replaced by a new element, whose initial opinion and affin- 
ity are randomly fished, with uniform probability, from 
the respective intervals [0,1] and [0,Q! max ]. Further, we 
introduce p — ^ to characterize the departure density, 
a crucial quantity that will play the role of the control 
parameter in our subsequent developments. As a final re- 
mark, it should be emphasized that no aging mechanisms 
are introduced: agents are mature enough to experience 
peer to peer encounters from the time they enter the sys- 
tem. 



3 Results 

Numerical simulations arc performed and the evolution of 
the system monitored. Qualitatively, the system shows the 
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typical critical behavior and asymptotic clusters of opin- 
ion, as observed in its original formulation [4j. However, 
peculiar distinctions are found, some of those being ad- 
dressed in the following discussion. 



* o 
.0+ 



0.01 
P 



First, an apparently smooth transition is also observed '* 

within this novel formulation, which divides the mono- % 

\ 

Ska 

and multi-clustered phases. Interestingly, the transition 10 «| 

10"' 10 10 1 

a/a 

point is now sensitive to the departure density p. To fur- 
ther elucidate this point, we draw in Figure 1 the average Fig- 2 . Main plot: Average number of clusters as function of 
number of observed clusters versus a rescaled temperature, the rescaled quantity f- for different values of the density p. 
A clear transition towards an ordered (single-clustered) Diamonds, stars, crosses, x-marks, circles and solid diamonds 
phase is observed, as the temperature increases. The pa- represent, respectively, simulations run at p = 0, p = 0.0025, 
rameter a c in Figure [2] plays the role of an effective tern- P = 0.0033, p = 0.005, p = 0.01 and p = 0.05. In all sim- 
perature, and it is numerically adjusted to make distinct ulations, here and after reported, unless otherwise specified, 

curves collapse onto the same profile, which hence holds Qc = AOc = °- 5 ' °° and a °'i are random variables uni- 

, £ ,i , c mi. ■ j. rw roi i formly distributed in the intervals [0,11 and [0,a maa .l - being 

irrespectively of the value of p. ihc inset of b igure[2| shows J L ' 1 1 ' 1 6 

that there is a linear correlation between a c and p. The 



otmax = 0.5 - respectively. Inset : a c as a function of p. The 
open circles represent the values calculated from the transi- 
tions shown in the main plot. The dotted lines represent the 
best linear fit : a c = 0.3p + a co , with a co = 5.5 • 10~ 4 . No- 
tice that <T C = <Jco is eventually recovered in the closed-system 
setting, which in turn corresponds to p — 0. 



larger the departure density p, the larger the effective tem- 
perature a c . In other words, when p is made to increase 
(i.e. the system is experiencing the effect of a more pro- 
nounced external perturbation), one needs to augment the 
degree of mixing, here controlled by the social temperature 

a, if a convergence to the final mono-cluster is sought. The , , „ n 

' b ° agents converge to a single cluster. When p is set to zero, 

death/birth process here modeled is in fact acting against , , n , , 1 n r ± 1 • 1 n 

1 r t> t> the final sfiared opinion is 0.5 to which all agents even- 

the thermal contribution, which brings into contact oth- t^. 1™ T , c , 

' ° tually agree, see figure |3p,. In other words, the final dis- 

erwise socially distant individuals. While this latter effect , . u t-v- j u .,1 ,1 1 • j • 

J tribution is a Dirac delta, with tne peak positioned m 

enhances the chances of convergence, the former favors the ^ n - ^ , c ^ n c n j- , ,• c 

b ' U = 0.5. Conversely, for p > 0, the final distribution of 

opposite tendency to spread. . . , , , , , , 

opinions presents a clear spreading around tne most prob- 

To further elucidate the role of the external pertur- able value, still found to be 0.5. This scenario is clearly 

bation, we shall refer to the dynamical regime where the depicted in Figure [3p. 



6 



Timoteo Carletti et al.: Opinion Dynamics in an Open Community 



0.6 



0.2 



0.5 
Opinion 



1 _ 2 

Time 



x 10 



0.8 



0.6 




1 2 

Time 



x 10 



Fig. 3. Opinion evolution versus time, the latter being quantified through the number of iterations. The black stars represent 
the mean opinion and the error bar the standard deviation of the opinion distribution. In the insets, the histogram of asymptotic 
distribution of agents' opinion. In panel a) the system is closed (i.e. p = 0), in b) p = 5 ■ 1CP 3 . 



The associated standard deviation v is deduced, from 
a series of simulations, and shown to depend on the se- 
lected value of p. The result of the analysis is reported 
in Figure EJ where the calculated value of v (symbols) is 
plotted versus the departure density amount p. 

For small values of the control parameter p, the stan- 
dard deviation v of the cluster scales proportionally to 
y/p, Figure [3] Numerics indicates that the proportional- 
ity coefficient gets smaller, as a c grows. In the opposite 
limit, namely for large values of the density p, the stan- 
dard deviation v rapidly saturates to a asymptotic value, 
v c . The latter is universal, meaning that it neither scales 
with a c , nor it does with the social temperature a. Our 
best numerical estimates returns v c = 0.28. 

The solid lines in Figure [5] represent the function : 

M 



which straightforwardly follows from a simple analytical 
argument, developed hereafter. In the above expression T c , 
stands for the time of convergence of the opinion cluster. 
In [6], working within the Dcffuant's scheme [3], i.e. in the 
closed community case without affinity, authors showed 
that the time needed to form a coherent assembly from a 
sequence of binary encounters scales as NlogN, where N 
refers to the population size. In 0], it was also proven that 
the affinity slows down the convergence rate and a correc- 
tion for T c that accounts for the role of a c was worked 
out. More precisely, it was found that : 



12^1-^(^)2]' 



(6) 



T c (a c ) 



N\og{N) 



21og(2) (a*-Oo) ' 



(7) 
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approximated form: 




Fig. 4. Main: The standard deviation of the final mono-cluster 
v as a function of the departure density p. Each point results 
from averaging out 100 independent runs. Symbols refer to 
numerical simulation performed with different values of a c , M 
and a (stars: a c — 0.5, M = 1, a = 0.07, diamonds: a c — 0.5, 
M = 2, a = 0.07, squares: a c = 0.5, M = 6, a = 0.07, dots: 
a c = 0.5, M = 2,a = 0.25, down - triangles: a c = 0.5, M = 2, 
cT — 1, open circles: a c = 0.3, M — 2, a — 0.28, up - triangles: 
q c = 0.3, M = 6, a = 0.28 and crosses: a c = 0, M = 2, 
cr = 0.28. The solids line refers to the theoretical prediction ((6]). 
Inset: v versus yfp, for p £ [0, 0.5]. This zoomed view confirms 
the correctness of the scaling dependence derived in ([8]). 



where a* ~ 0.64[J. Before turning to discuss the analytical 
derivation of Eq. ([6|), we wish to test its predictive ade- 
quacy with reference to the two limiting cases outlined 
above. Indeed, for p « 1, Eq. (J6]) can be cast in the 



3 a, is an effective value deduced via numerical fit, as dis- 
cussed in [4]. A fully analytical derivation of T c will be pre- 
sented in [7]. 



T c (a c 



<M 



(8) 



24iV V T 

which presents the same dependence of v versus y/p, as 
observed in the numerical experiments. Moreover, the co- 
efficient 7 t is predicted to decay when increasing the cut- 
off in affinity a c , just as expected from the numerics. For 
p » 1, Eq. ([6]) implies: 

v = ^^0.28, (9) 

in excellent agreement with our numerical findings. 

To derive Eq. © let us suppose that at time t the 
death/birth process takes place and the system experience 
an injection of new individuals. Label with v t the stan- 
dard deviation of the agents opinion distribution / (O), 
at time t. It is reasonable to assume that /*(0) is cen- 
tered around 1/2. After T interactions between agents, 
when the next perturbation will occur (M agents are ran- 
domly removed from the community and replaced by M 
new actors with random opinion and affinity scores) the 
distribution has been already modified, because of the un- 
derlying dynamical mechanism specified through Eqs. ([3]). 
More concretely, the opinions slightly converge around the 
peak value 1/2, an effect that certainly translates into a re- 
duction of the associated standard deviation. To provide a 
quantitative estimate of such phenomenon, we recall that 
in the relevant (O, t) plan, the convergence process fills an 
ideal triangular pattern, whose hight measures T c . This 
topological observation enables us to put forward the fol- 
lowing linear ansatz: 

„„„., T r - T 



J t+T 



T 



(10) 



8 Timoteo Carletti et al.: Opinion Dynamics in an Open Community 

where labels the standard deviation just before the state, being characterized by a single cluster, centered in 

insertion of the next pool of incoming agents. Recalling O — 0.5 and with standard deviation given by Eq. [6j It is 

that the newly inserted elements are uniformly distributed, also observed that the time needed by the system to com- 

and labeling with G(-, •) the opinion distribution (the two plete the transition T conv depends on the value of p and 

entries referring respectively to mean and the standard the critical affinity a c , see Figure [5] A simple theoretical 

deviation), the updated variance is: argument enables us to quantitatively explain these find- 



M r 1 / ] \ 2 ings. The initial distribution of opinions is modified after 

the first death/birth event as: 



N - M ' 1 g(1^\Io-\\\o 



^» m _- + ^ (0 _o.>. (11) 

M N-M (T r — T x iv iv 



127V N V T 







The first term refers to the freshly injected actors, while 

The asymptotic stationary solution correspond to t ^ e sec0 nd stands for the remaining Delta-distributed in- 

v t +T =v t ,& condition that immediately leads to Eq. © dividuals. Hence, the mean opinion value reads: 

when plugged into (TTT1) . The above analysis also suggests _ r 1 , M N — M 

Oi = I f\0)OdO = — + ^^0 . (12) 

that the final fate of the system is not affected by the 
time when the perturbation is first applied, t i n . This con- 
clusion is also confirmed by direct numerical inspection: 
The asymptotic value of the standard deviation v does 
not depend on i,- n , but solely on p. Even in the extreme 
condition, when the death/birth perturbation is switched 
on after the agents have already collapsed to the mean 
opinion 0.5, one observes that, after a transient, the clus- 
ter spreads and the measured value of v agrees with the 

theoretical prediction ([6]). 

as: 

Aiming at further characterizing the system dynam- 

d*=*L + H^±di. (13) 
ics, we also studied the case where, initially, agents share 2N N 

the same belief O . The initial distribution of opinions is Aftcr 71 death/birth iterations, the opinion mean value 
therefore a Dirac delta f°(0) = 6(0 - D ). Such con- reads : 

- _ M_ N-M ^ 

dition is a stationary solution, for any given when ^" — 27V N ™~ 1 ' ' > 

the death/birth process is inactivated. Conversely, when From Eq. (TT?)) one easily gets that the asymptotic equilib- 
the death/birth applies, the system evolves toward a final rium is reached for Ooo = 0.5, as found in our numerical 



2N N 

We can suppose that between the occurrence of two 
consecutive perturbations (separated by T iterations), the 
group average opinion does not significantly change. No- 
tice that the probability of interaction of a newborn agent 
with another belonging to the main group is in fact pro- 
portional to M/N. Moreover several consecutive encoun- 
ters of this type are necessary to induce a macroscopic 
change of the averaged opinion. Under this hypothesis the 
next death/birth event makes the average opinion change 
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experiments; in fact the following relation is straightfor- 
wardly obtained: 



Or, 



N -M 
N 



' \A M (N-M\ 



1=0 



N 



being O the initial common believe. By setting a = 
and (3 = Mr, the solution of Eq. (fT5")) reads: 



1 - a n 

O n = (3- + a n O a , 

1 — a 



N-M 
N 



(16) 



whose asymptotic solution is given by n — > Ooo = 0.5. 

Expression (|16|) reproduces quite well the dynamics of 
the cluster mean, as seen in the simulations. The adequacy 
of (flo) is in fact clearly demonstrated in Figure [5^. Let 
us define the convergence time T conv as the number of 
iterations needed to bring the average opinion e close to its 
asymptotic value 1/2. Solving Eq. (|16p for n and recalling 



that T,., 



nv = nT yield: 

Tconv = T log Q 



TT( e - 1 ) 

77 1 — a 



(17) 



\Oo 2 

The above estimate is in excellent agreement with the nu- 
merical results reported in Figure Et>. 
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Fig. 5. a) The opinion as a function of time for the case where 
all the agents are initialized with the same opinion O = 0.1. 
The solid (red) line represent the numerical data while the 
segmented (blue) line represents the theoretical prediction from 
Eq. (|16p . b) The transition time T conv as a function of the death 
density p. Symbols stand for the results of the simulations and 
the solid line represents the prediction 1)170 . Here iV = 100, 
a c = AO c = 0.5, e = 0.001 and a = 0.07. 



4 Conclusions 

In this paper we have discussed the process of opinion 
making in an open group of interacting subjects. The 
model postulates the coupled dynamical evolution of both 
individuals' opinion and mutual affinity, according to the 
rules formulated in [3]. At variance with respect to the 
toy — model [4j, the system is now open to contact with 
an external reservoir of potentially interacting candidates. 
Every T iterations the M agents are instantaneously re- 
placed by newborn actors, whose opinion and affinity scores 



are randomly generated according to a pre-assigned (here 
uniform) probability distribution. The ratio p = M/T, 
here termed departure density plays the role of a con- 
trol parameter. The occurrence of a phase transition is 
found which separates between two macroscopically dif- 
ferent regimes: For large values of the so-called social 
temperature the system collapses to a single cluster in 
opinion, while in the opposite regime a fragmented phase 
is detected. The role of p is elucidated and shown to en- 
ter in the critical threshold as a linear contribution. Two 
phenomena are then addressed, with reference to the sin- 



10 Timoteo Carletti et al.: Opinion Dynamics in an Open Community 

gle clustered phase. On the one side, the external pertur- 
bation, here being hypothesized to mimic a death/birth 
process, induces a spreading of the final cluster. The as- 
sociated variance is numerically shown to depend on the 
density amount p, the functional dependence being also 
analytically explained. On the other, we also show that 
the birth/death events imposed at a constant pace can 
produce the progressive migration of a cluster, initially 
localized around a given opinion value. A theoretical argu- 
ment is also developed to clarify this finding. As a general 
comment, we should emphasize that the effect of opening 
up the system to external influences changes dramatically 
its intrinsic dynamics revealing peculiar, potentially inter- 
esting, features which deserves to be further explored. 
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